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ABSTRACT 

The millimeter-wave polarization of the protoplanetary disk around HL Tau has been interpreted 
as the emission from elongated dust grains aligned with the magnetic field in the disk. However, the 
self-scattering of thermal dust emission may also explain the observed millimeter-wave polarization. 
In this paper, we report a modeling of the millimeter-wave polarization of the HL Tau disk with the 
self-polarization. Dust grains are assumed to be spherical and to have a power-law size distribution. 
We change the maximum grain size with a fixed dust composition in a fixed disk model to find the 
grain size to reproduce the observed signature. We find that the direction of the polarization vectors 
and the polarization degree can be explained with the self-scattering. Moreover, the polarization 
degree can be explained only if the maximum grain size is ~ 150 /im. The obtained grain size from 
the polarization is different from a size of millimeter or larger that has been previously expected 
from the spectral index of the dust opacity coefficient if the emission is optically thin. We discuss 
that porous dust aggregates may solve the inconsistency of the maximum grain size between the two 
constraints. 

Subject headings: dust, polarization, protoplanetary disks 


1. INTRODUCTION 

Protoplanetary disks are believed to be the birth¬ 
place of planets. Due to the high spatial dust den¬ 
sity, submicron-sized dust grains coagulate to form larger 
bodies, and ultimately form planets in the disks (e.g. 


Weidenschilling 


1977 Hayashi 19811. Thus, observa¬ 


tional constraints on the grain size in protoplanetary 
disks are essential to directly investigate the ongoing 
planet formation. 

The grain size in protoplanetary disks has been con- 


ficient at millimeter wavelengths (e.g., 

Beckwith et al. 

1990 

Beckwith & Sargent 1991 

Miya 

kb & Makagawa 

1993 

Andrews & Williams 2005; 

Iselia et al. 2009; (Ricci 

et al. 

2010b |a van der Maref et al. 2013). 4'he opacity 


has been shown to be as low as 0 — 1, which indicates 
that the grain size in protopl anetary disks is of the order 


of millimeter, or larger (e.g., Draine 2006). In addition, 
even at some of the envelopes of protostars, the opacity 
index is indicated to be as low as the later stage of pro¬ 
toplanetary disks, which suggests that the dust grains 
are grown to millimeter i n size even in the ear ly phase of 
the protoplanetary disks (Miotello et ah][2014) although 


the modeling contains large uncertainties. However, the 
low opacity inde x can also be expl ained with optically 


protoplanetary disks. 


kataoka@uni-heidelberg.de 


Recently, an independent method to constrain the 
grain size in p rotoplanetary disks has been proposed. 
|Katao k a et al.j2015| proposed that millimeter-wave emis- 
sion is partially polarized due to the self-polarization if 
the following two condition are satisfied. The first con¬ 
dition is that dust grains have the size comparable to 
the wavelengths to have the scattering efficiency large 
enough to produce the scattered emission. The second 
is that the thermal dust emission has anisotropic distri¬ 
butions so that the thermal dust emission is scattered 
by dust grains themselves and to show the residual po¬ 
larization in total due to the anisotropic distribution of 
the incoming fluxes. Therefore, if the polarization due to 
the self-scattering of dust thermal emission is detected, 
it would be the evidence of the existence of dust grains 
which have the comparable size to the wavelength. In 
this paper, we report the application of this method to 
the protoplanetary disk around HL Tau. 


HL Tau is a p rotostar 140 pc away from the sun (Re¬ 


bull et al. 2004). This young star is surrounded by an 
envelope and producing jets and outflo ws, which indi¬ 


cates a highly active star formation (e.g., Hayashi et al. 
1993). The millimeter emission from the envelope and 


the disk ha s been intensively investigated with inter¬ 
ferometers ( Beckwith et al. 19901 Mundy et al. 1996 


thick disks 
grains (e.g., 
cal composit 

e.g., Ricci et al. 

2012 

irregularly shaped 
or different chemi- 

Min et al. 2005 

2007) 

ion of the grams 

Al¥) 

Pollack et al. 1994 

Mennella et al. 1998 Jones 2012). 

Therefore, there are 


Looney et af.||2U00||Wifner et af.||1996[ Greaves et al. 

2008[ Kwon et al. 2011). Furthermore, Atacama Large 


Millimeter/submilhmeter Array (ALMA) has revealed 
that the circumstellar disk around HL Tau has multi¬ 
ple rings at sub-millimeter wa veleng ths, which mig ht be 
a sign ature of multiple planets (ALMA Partnership et al. 


20151. In addition, spatially resolved polarized emission 
of the disk at sub-nun wavelengths has been detected 
with Combined Array for Millimeter-wave Astronomy 
(CARMA) and Submillimeter Array (SMA). The net po- 
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larization degree is 0.89% in average with CARMA, and 
0.86 ± 0.4% with SMA. The polarization vectors are di¬ 
rected from north-east to south-west. In this paper, we 
construct a disk model which reasonably reproduces the 
ALMA observations and perform radiative transfer cal¬ 
culations to investigate the polarization due to the self¬ 
scattering of thermal dust emission. 

The millimeter-wave polarization observations toward 
star-forming regions and protoplanetary disks have been 
interpreted as an indicator of the morphology of the mag- 
netic field (Girart et al.|2006| 2009 Hull et al.|2013 2014 


Rao et ah] 2014| |Stephens et al.||2014| |begura-(Jox et al. 

2015[). The elongated dust grains have a preferential di- 


rection perpendicular to the magnetic field, which results 
in the polarization of the thermal dust emission (Davis 


& Greenstein 1951). In addition, radiative torq ues help 
dust grains to be aligned with the magnetic field (Draine 


fe Weingartner|1997l|Cho fe Lazarian|2005[ [2007| | Hoang 


& Lazarian|[2008 2009b|a|). 

This paper does not intend to exclude the possibility 
that the polarized emission of HL Tau disk is caused 
by the grain alignment with the disk magnetic field but 
investigates the self-polarization as an alternative expla¬ 
nation for the polarization observed with CARMA and 
SMA. We also discuss possible methods to distinguish 
the mechanisms between the grain alignment and the 
self-polarization. 

Note that , during the preparation of the manuscript, 
|Yang et ah1 [2015 independently found that the interpre¬ 
tation of the millimeter-wave polarization of HL Tau due 
to dust scattering. 


2. METHOD 

We construct an axisymmetric dust disk model with 
several gaps with a smooth temperature distribution 
which reasonably repro duces the sub- mm continuu m im¬ 
age at A = 1.3 mm (ALMA Partnership et al. 2015). 
This model is not a unique solution of the density and 
temperature of the HL Tau disk. However, this model 
would be enough to investigate the grain size with po¬ 
larization signature as long as the millimeter-wave con¬ 
tinuum is reasonably reproduced. The polarization frac¬ 
tion is determined by the combination of the grain size 
and the anisotropy of the thermal emission. The differ¬ 
ence in the temperature and density distributions affects 
the anisotropy of the radiation field, but the dependence 
of the polarization on the grain size is t houg ht to be 
stronger. We discuss this point in Section |4.4| but the 
detailed parameter studies of the density and tempera- 
ture remain to be a future work. 

First, we fix the temperature profile with the smooth 
power-law distribution as 


T = T 0 



( 1 ) 


*£ 


^1 — exp 


( I f R - r d ,i 

\ 2 \ W d,i 



The parameters are summarized in Tables [T] and [2j 



Fig. 1.— The data points show the observed intensity profile of 
Band 6 of ALMA along the major axis in brightness temperature. 
The solids line shows the brightness temperature of the disk model 
adopted in this paper. 


Figure [I] shows the brightness temperature derived 
from the Band 6 data on the major axis illustrated with 
the model brightness temperature, which is calculated as 

I V {R) = (1 - exp(— t(.R)) B v {T{R)). (3) 


Then, we model the disk with dust surface density of 
Ed — T(l?)/ftabs (^max) A — 1.3 mm), where ^abs(^max) A) 
is the absorption opacity including a size distribution of 
dust grains. 

In calculating the opacity, the dust grains are assumed 
to be spherical and have a p ower -l aw size distr ibution 
with a power of q = —3.5 (Mathis et al. 1977) with 
the maximum grain size a max . We take this maximum 
grain size as the representative grain size in the follow¬ 
ing discussion. The opacity is calculated with Mie the¬ 
ory. The composition is assumed to be the mixture of 


silicate, organics, and water ice (Kataoka et al. 2014 
Pollack et al.|1994). We use the retractive index of astro¬ 


nomical silicate (We ingartner & Draine 2001), organics 
(Pollack et al. 1994|), and water ice (| Warren 


1984) and 

calculate the mixture of them with the effective medium 


theory with the Maxwell-Garnett rule (e. g., Bohren & 


Huffman 1983 Miyake & Nakagawa 1993). A different 


abundance maylead to the different absolute value of 
polarization degree, which should be investigated in fu¬ 
ture studies. The adopted value for the fiducial run is 
«abs(a m ax = 150 /rm, A = 1.3 mm) = 0.24 cm 2 g -1 . In 
the fiducial case, therefore, the dust surface density Ed 
has a profile of 


Then, we add the additional gaps on the density pro¬ 
file to reproduce the observation along the major axis of 
the disk observed with the continuum observed with the 
ALMA Band 6. 


T = T 0 



R 


Re 




( 2 ) 


E d = 8.3 g cm 2 


R 


1 AU 



*£ 


^1 — exp 


( _ 1 / R - r d ,i 

y 2 \ w dji 



(4) 


















































































3 


TABLE 1 

The disk parameters 


parameters 

values 

To 

280 K 


0.3 

TO 

2.0 

p 

0.3 

^exp 

4 

Rexp 

120 AU 


TABLE 2 

The gap parameters 


gap number 

r d,i 

Wd,i 

fd,i 

1 

12.5 

3.0 

1.25 

2 

32.0 

4.0 

1.0 

3 

42.0 

2.0 

4.0 

4 

55.0 

4.0 

3.0 

5 

65.0 

4.0 

1.4 

6 

77.0 

6.0 

1.5 

7 

95.0 

5.0 

2.0 


We perform the radiative transfer simulations with a 
public code RADMC-3D to obtain the dust continuum, 
the polarized intensity, and the polarization degree. To 
obtain the vertical density distribution, we assume the 
Gaussian density distribution with a dust scale height 
hd such that pd = Sd/(V27r^d) exp(— z 2 /2h\). To repro¬ 
duce the geometrically flat disk observed with ALMA, 
we set hd = h g //settle > where h g is the gas pressure scale 
height and / se ttie = 10 (see Appendix). Here, we do not 
use the thermal Monte Carlo simulations to determine 
the temperature but use the simple power-law temper¬ 
ature model described above. We assume that the dis¬ 
tance to HL Tau is 140 pc, so 1 arcsec corresponds to 
140 AU in the figures. The inclination is assumed to be 
40°. 

Note again that this modeling is not a unique solution 
to reproduce the emission of the HL Tau disk. However, 
the main goal of this paper is to constrain the grain size 
from the polarization observations. As we discuss in the 
following sections, the polarization degree is mainly de¬ 
termined by the combination of the observed wavelengths 
and the grain size, which does not so much depend on the 
detailed modeling of temperature and the surface density. 


3. RESULTS 

We show the results in the case that the maximum 
grain size is a max = 150 pm as a fiducial case. Figure [2] 
(a) shows the dust continuum, where we confirm that the 
continuum image well reproduces the multiple-ring struc¬ 
ture observed with ALMA. Figure [2] (b) shows that the 
polarized intensity overlaid with the polarization vectors. 
Figurej2](c) shows that the polarized degree overlaid with 
the polarization vectors. 

The basic feature of the polarized intensity and the 
polarization vectors can be explained with the self- 
polarizati on with the anisotro py of the thermal dust 
emission (Kataoka et al. 2015). The thermal emission 
from dust grams are originally unpolarized. The unpo¬ 
larized light is scattered by other grains because of the 
high scattering opacity. Due to the ring-like structure of 
the emission and the inclination of the disk, the distri¬ 


bution of the incoming fl ux to t h e scatterin g dust grains 
has anisotropies (see also Yang et al(| |2015 1. As a result, 
the total flux has a residual polarization corresponding 
to the anisotropic radiation field. Note that the scat¬ 
tered emission comes mainly from midplane because the 
disk is optically thin or marginally thick in the vertical 
direction (see Equation (J2|). 

The polarized intensity is centrally concentrated. This 
infers that if we detect the polarization with a marginal 
sensitivity, we can detect the central part of the disk. 
This is consistent with the results of polari zation obser¬ 
vation with CARMA ( Stephens et al. 12014 ), which show 
the centrally concentrated polarized intensity. In addi¬ 
tion, the polarization vectors are directed from top left to 
bottom right (north-east to south-west). This is also con- 
sistent with the polari zation observation with CARMA 
(JStephens et al.| 2014|. Figure [2| (c) shows the polar¬ 
ization degree overlaid with polarization vectors. The 
polarization degree is two times higher in the gap re¬ 
gions than that in ring regions. At present, the polariza¬ 
tion image of CARMA does not have a spatial resolution 
high enough to resolve the rings. Future observations of 
ALMA with sp atial resolu tion as hig h as t he lon g base¬ 
line campaign (ALMA Partnership et al. 2015) would 
reveal these structure even m the polarization. 

The reinterpretation of the polarization puts a strong 
constraint on the size of dust grains. Figure [3] shows 
the theoretical model of the net polarization degree at 
the wavelength of A = 1.3 mm as a function of the maxi¬ 
mum grain size a max . The top of Figure [3] shows the total 
polarization degree for several radiative transfer calcula¬ 
tions at A = 1.3 mm as indicated black dots overlaid 
with the theoretical model. The theoretical model is cal¬ 
culated as the product of pola rization degree at 90° scat¬ 
tering P go and the albedo oj (Kataoka et al. 2015). The 
expected polarization degree is 


P = CP 90 uj, (5) 

where C = 2.0 % calibrated with the radiative transfer 
calculations in this paper. We also includes the u ncer- 
tainty of the model of ± 50 % (see Appendix of|Kataoka 
|et aL| |2015| for the error of the polarization degree in 
different radiative transfer calculation codes). The ob¬ 
served polarization degree with CARMA is also shown. 
This figure clearly shows that the observed polarization 
degree of 0.89 % at A = 1.3 mm observed with CARMA 
can only be explained with the grain size in the range of 
70 /im < a max < 350 /im. If the grain size is much lower, 
the scattering opacity is too low to scatter the thermal 
dust emission. If the grain size is much higher, the scat¬ 
tering is forwardly peaked and thus no polarization can 
be expected. 

We put further constraints on the grain size with the 
results of SMA. The polarization degree is a function 
of the combination of the observed wavelength and the 
grain size. Therefore, the results of the different wave¬ 
lengths put further constraints. The bottom of Figure 
[3] shows the expected polarization degree as a function 
of observed wavelengths with changing the maximum 
grain size a max . As shown in the figure, the result of 
Umax = 350 /iin explains the CARMA observations at 
A = 1.3 mm but underestimates SMA observations at 
A = 0.89 mm. Therefore, we can rule out the possibility 
of a max = 350 /im to explain the polarization. In this 
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Fig. 2.— (left) The intensity map of the radiative transfer calculations in the unit of [Jy/arcsec 2 ]. (center) The polarized intensity 
map in the unit of [mJy/arcsec 2 ], overlaid with the polarization vectors, (right) The map of the polarization degree in the unit of %, 
overlaid with the polarization vectors which is the same as the middle panel. 



way, we put further constraints on the grain size to be 
in the range of 70 /rm < a max < 270 p,m to explain the 
polarization degree of the both CARMA and SMA ob¬ 
servations. We take a max = 150 /rm as a representative 
value and continue the discussion. 

4. DISCUSSION 

4.1. Grain size constraints with opacity index 

The grain size has been constrained with the spectral 
index at sub-mm wavelengths. If the emission is optically 
thin, the index has the information of the dust opac¬ 
ity index at sub-mm wavelengths. The opacity index is 
typically equal to or less than 1 in protoplanetary disks, 
which can be explained wit h millimeter-sized grains (e.g., 
Beckwith & Sargent 1991). In the case of HL Tau, the 
dust opacity ind ex is m the range from 0.3 to Oil for 
the bright rings (ALMA Partnership et al. 2015 Kwon 
et al. 2011). Therefore, the maximum gram size of 150 
pm, which is obtained in this paper, is not consistent the 
constraints on the grain size with the opacity index if the 
emission is optically thin. 

Figure [4] shows the opacity index (3 in the case of 
q = —3.5 and q = —2.5 in the case of the dust model 
adopted in this paper. We change the power-law index 
because there is no solution if the power of the grain 
size distribution is q = —3.5. To explain the observed 
opacity index, the grain size should be in the range of 
1 cm < a max < 5 cm if the power of the grain size distri¬ 
bution is q = —2.5. The opacity index strongly depends 
on the composition of the dust grains. For example, dust 
grains composed of silicate and carbonaceous materials 
can produce the lower value of beta than the ice-included 
grains; the grains with mixture of silicate and carbona¬ 
ceous material c an reach 0 = 1 ev en with a max = 500 pm 
(see Figure 4 of Testi et al.|2014 for example). Although 
there are uncertainties m compositions, the maximum 


grain size expected from the interpretation of the spec¬ 
tral index is significantly larger than that expected from 
the polarization. 

4.2. Porous dust aggregates as a possible solution 

Here, we discuss the porous dust aggregates as a pos¬ 
sible solution to solve the inconsistency of the grain size 
between the two constraints. We have constrained the 
size of maximum dust grains under the assumption of 
spherical grains. The upper limit of the maximum grain 
size is determined by the fact that the large spherical 
dust grains do not show the polarization due to scat¬ 
tering but forwardly scatter the light. However, if we 
consider porous dust aggregates, the properties of their 
scattering and resultant po larizati on reflect the p roper¬ 


ties of constituen t particles (Min et al. 2015 Tazaki et al. 
submitted to ApJ), which are believed to have the size of 


(sub-)micron size. Therefore, if we consider dust aggre¬ 
gates which have a even larger size than the constraint 
from the millimeter-wave polarization, they may explain 
the polarization. This infers that the highly porous and 
massive aggregates may have both of the low opacity in¬ 
dex and the high efficiency of scattering to produce the 
polarization. 

To understand what kind of dust aggregates are re¬ 
quired to explain both of the polarization and the spec¬ 
tral index, here we discuss the possible constraints on 
the dust aggregates from the spectral index assuming 
that the emission is optically thin. Dust g rains coag¬ 
u late to form porous dust aggregates (e.g., Ossenkopf 


19931 The filling factor can be even as small as 10 4 in 


disks (Kataoka et al. [2013 ) although how porous is the 
dust aggregates is still controvercial. However, we can 
constrain the product of the aggregate radius a and the 
filling factor / because the abs orption opacity is the same 
if the product af is the same (Kataoka et al. ]M4l ). Fig- 
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Fig. 3.— (top) The net polar ization degrees of 0.90% (average) and 0.72% (median) with CARMA observations are shown as red and 
green lines ||Stephens et al.|2014|). The red shaded region is the analytical model of the net polarization degree. The black points represent 
the net polarization degree ot each radiative transfer calculations. The analytical model has a good agreement with radiative transfer 
calculations. From this figure, the grain size is constrained to be 70 /im < a max < 350 //m. (bottom) The net polariz ation degrees of 
0 .90% (average) and 0.72% (median) with CARMA observations and 0.84% ± 0.4% with SMA observations are shown ||Stephens et ah 
(2014). Expected polarization degree as a function of observed wavelengths. Each thick line represents the several cases oi the gram size, 
which are a max = 50, 70,150, 270, 350/un. The lines include the error of ±50%. 


ure [4] also shows the opacity index as a function of the 
product of grain/aggregate size a and the filling factor 
/. In the case of q = —2.5, the grain size to explain 
the observations is 5 cm < a max f <21 cm in the case of 
the fluffy aggregates (/ < 0.1). To explain the observed 
index of the absorption opacity, the dust aggregates in 
HL Tau should have a product of af in the range of 
5 cm < a max / < 21 cm in the adopted dust model. This 
means that, for example, if the filling factor is as small 
as / = 10~ 4 (Kataoka et ah][2013), the aggregate radius 
is ~ 1 km. 

The size distribution does not significantly affect the 
polarization degree. Figure [5] shows that the expected 
polarization degree in different cases of the power of the 
size distribution. We do not see a significant effects on 
the polarization. Therefore, we conclude that the the 
grain size constrained with the polarization is more ro¬ 
bust than that with the opacity index. 


4.3. Stokes number 


The dynamics of dust grains or aggregates are deter¬ 
mined by the normalized stopping time due to the fric¬ 
tion between the gas and dust, which is the Stokes num¬ 
ber. The great advantage of knowing the product af is 
that we can determine the Stokes number if the aggre- 
gate radius is smaller than the mean free path of the gas 
(Kataoka et al. 2014). If we obtain the Stokes number, 
we can discuss the coupling efficiency between the dust 
and gas. We will discuss the Stokes number based on 
the discussion above. Here, we consider two dust grain 
size: the 150 /im-sized compact grains, which explain the 
polarization observations but do not explain the spec¬ 
tral index, and the 1 km-sized fluffy aggregates with 
/ = 10~ 4 , which explain the opacity index and which 
may explain the polarization properties although it is 
highly uncertain. The Stokes number is written in the 
form of St ~ p mat a/E g . With the adopted model, the 
Stokes number is estimated to be St ~ 10 -4 — 10 -3 in 
the case of the 150 //m spherical grains. If the dust ag¬ 
gregates are fluffy and they have af = 10 cm, St ~ 10” 2 . 
This means that if we assume compact grains, the possi- 
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Fig. 4.— (left) The opacity index between ALMA Band 6 (1.3 mm) and Band 7 (0.89 mm) as a function of the product of the maximum 
grain size a max and the filling factor. The power of the size distribution is taken to be q = —3.5. The observed range of the opacity index 
is indicated as shaded region, (right) The power of the size distribution is taken to be q = —2.5. 



Fig. 5. — The expected polarization of the circumstellar disk of 
HL Tau with different grain size distribution. The red and green 
solid lines represent CARMA observations (see Figure pjf. We take 
a power-law grain size distribution of n(a) oc a q and change q as 
q = —3.5, —3.0, —2.5, —2.0. 


gt al. 

2009), or a 

|ierro 

et al. 

2015) 


tive velocity compared with the compact grains. Future 
observations on the gas density distribution at the rings 
will reveal them. 

In addition, the geometrical thickness of the disk also 
hints the grain size. As shown in the Appendix, the 
observed dust continuum image can no t be re pr oduce d 
without the vertical settling (see also Pinte et ah]|2015). 
The vertical thickness of the disk is considerably smaller 
than the thermal scale height of the gas, which indi¬ 
cates that the dust grains are settled toward the mid¬ 
plane. This requires at least the Stokes num ber i s la rger 


than the turbulent parameter a (e.g., Youdin & Lithwick 


2007). Therefore, the turbulent parameter a should at 
least be lower than 10 -4 — 10 -3 in the case of compact 
grains and be lower than 10 -2 in the case of fluffy dust 
aggregates from on the constraints on the Stokes number 
we discussed above. 


ble grains which can reproduce the polarization is tightly 
coupled to the gas. If we assume that fluffy dust aggre¬ 
gates, the coupling to the gas is relatively weak. 

In the case of the spherical grains, the Stokes number 
is so small that they are coupled to the gas. Therefore, 
if it is the case, the gas and dust have the same density 
distribution as the multiple-ring structure, which enco ur¬ 
ages t he sc enarios of gas clearance by plane ts (|Tamayo 


et al. 


Z0T5 


2015), magnetic instabilities o f the gas (j ltockjetah 


, or gravitational instability (Takahashi & Inutsuka 


20l4). Furthermore, the rings may be related to en¬ 
hance d grain growth at snow lines of several dust spices 
(Zhang et al. 2015), which should be consistent with 
the spherical gram size because the condensation growth 
leads grains to be larger spherical grains. 

On the other hand, if the sub-mm emission is com¬ 
ing from the fluffy aggregates, the Stokes number may 
be large enough to be decoupled from the gas. This en¬ 
courages the sce narios of trapping the dust at a radial 


pressure bump (Pinilla et al. 2012), at a vortex (Lyra 


4.4. Dependence on disk models 

In this paper, we have used only one disk model to 
constrain the grain size by the polarization signature. 
However, the dependence of the polarization fraction on 
the disk model is weak. As shown in Figure [3] the po¬ 
larization fraction is described as P = CP q ciQJ and P qquj 
is determined by the grain size (Kataok a et al.| |2015). A 
different disk model only changes the calibration parame¬ 
ter C, which is calibrated to be 2.0% in the adopted disk 
model in this paper. Therefore, unless the disk model 
can change the calibration parameter C in a factor of 
few, the constraint on the maximum grain size in this 
paper does not change. 


4.5. How to distinguish the polarization mechanisms? 

In this paper, we focus on the interpretation of the 
observed mm-wave polarization properties with the self¬ 
scattering mechanism. However, elongated dust grains 
aligned with the magnetic field can also explain the polar¬ 
ization. Here, we discuss the ways to distinguish the two 
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mechanisms. The most promising way to distinguish the 
polarization mechanisms is to perform the polarization 
observations at other wavelengths. Although the wave¬ 
length dependence of the polarization degree is strong in 
the case of the self-scattering as shown in Figure [3} the 
wavelength dependence in the case of the grain alignment 
is not so m uch strong although it is still uncertain (see 
Figure 2 of Andersson et al.||2015 for example). Further 
observations at wavelengths other than 0.87 mm and 1.3 
mm will give us a clue to distinguish the two mechanisms. 

5. CONCLUSIONS 

The protoplanetary disk around HL Tau s hows th e 
polarized emission at millimeter wavelengths (Tamura 
et al. 1995 Stephens et al. 2014). The polarized emis¬ 
sion has been interpreted as the thermal dust emission 
from elongated dust grains aligned with the magnetic 
field. However, the self-scattering of d ust g r ains m ay 
also explain the observed polarization (Kataoka et al. 
20151. Therefore, we have investigated whether the self- 


scattering of the thermal dust emission accounts for the 
observed millimeter-wave polarization of the protoplan¬ 
etary disk around HL Tau. We used a simple dust disk 
model which reasonably reproduces the millimeter-wave 
continu um observed w ith A LMA at 1 .3 millimeter wave¬ 
length (ALMA Partnershi p et alT] 2015). Dust grains are 
assumed to be spherical and have a power-law size dis¬ 
tribution with the power of q = —3.5. The maximum 
grain size a max is the parameter. We have performed 
radiative transfer calculations with the model described 
above with RADMC-3D to investigate the polarization 
properties at millimeter wavelengths. 

As a result, we successfully reproduced the polarization 
vectors and polariza tion degree of HL Tau observed with 
CARMA and SMA(Stephens et al. 2014). We changed 


the maximum grain size a max as a parameter and calcu¬ 
late the polarization degree to constrain the grain size. 
We found that the observed polarization degree can be 
reproduced only if the maximum grain size is in the range 
of 70 /mi < a max < 270 /xm. This is a strong constraint 
on the grain size in the protoplanetary disk around HL 
Tau. 

If the grain size is around 150 /im, it gives a constraint 
on the scenario of the trapping of dust grains through to 
the coupling efficiency between the gas and dust. Stokes 
number of these grains is estimated to be around 1CU 4 < 
St < 10~ 3 , which indicates the dust grains are almost 
coupled to the disk gas. 

We also discussed the possibility that the dust grains 
are porous in HL Tau disk. The Stokes number inferred 
from the spectral index is as large as St < 10~ 2 if the 
emission is optically thin and if the dust aggregates are 
highly porous. If this is the case, the dust aggregates are 
marginally decoupled from the gas. However, due to the 
lack of knowledge about the millimeter-wave polariza¬ 
tion of porous dust aggregates, we could not discuss the 
polarization degree. The further theoretical constraints 
on polarization properties of porous dust aggregates are 
required in future studies. 

We have also discussed the possible way to distinguish 
the mechanisms between the self-scattering and the grain 
alignment. One possible way is to perform the multi¬ 
wave observations because the polarization degree ex¬ 
pected from the self-scattering shows strong dependence 


on the observed wavelength although the wavelength de¬ 
pendence in the case of magnetic field alignment is not 
so much strong. Therefore, further observations at the 
wavelengths different from 1 mm is required in future 
studies. 


This paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2011.0.00015.SV. ALMA is a part¬ 
nership of ESO (representing its member states); NSF 
(USA); and NINS (Japan), together with NRC (Canada), 
NSC, and ASIAA (Taiwan); in cooperation with the Re¬ 
public of Chile. The Joint ALMA Observatory is op¬ 
erated by ESO, AUI/NRAO, and NAOJ. This work is 
supported by MEXT KAKENHI No. 23103004 and by 
JSPS KAKENHI No. 15K17606 and 26800106 


APPENDIX 

VERTICAL SETTLING OF DUST GRAINS 


We also perform the radiative transfer calculations 
without the dust settling; / se ttie = 1- In this case, the 
dust grains are well mixed with the gas. The other pa¬ 
rameters are set to be the same as the fiducial run and the 
maximum grain size is set to be 150 /xm. Figure [6] shows 
the intensity distribution. The gaps are not clearly seen 
in this picture. This is because the scale height of dust 
grains are too high to reproduce th e geometrically t hin 
disk observed with ALMA (see also Pinte et al.]|20151. 
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Fig. 6.— The intensity map in the case without the vertical 
settling of dust grains. 


DEPENDENCE ON DISK MODELS 

In this paper, we have used only one disk model, which 
reasonably reproduce the intensity distribution of ALMA 
observation at A = 1.3 mm. Although the adopted disk 
model is not a unique solution, the results of the con¬ 
straints on grain size would not so much depend on 
the disk model. This is because the different model of 
the density and temperature distribution gives different 
anisotropic radiation field but it only changes the abso¬ 
lute value of the polarization degree but does not change 
the relative dependence of the polarization degree on the 
grain size. To demonstrate this, we perform a radiative 
transfer calculation of another disk model in this section. 
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Fig. 7.— This figure is the same as Figure [l]but add the profile 
of the shallower temperature model adopted in the Appendix. 


Fig. 8 .— The same figure as Figure [3] but for the shallower 
temperature slope disk model. The analytical model has the form 
of CPqqui where C = 1.3%. 

We change the power-law index of temperature from 
qt = 0.3 to q t = 0.5. In addition, to roughly fit the 
overall radial distribution of intensity, we set p = —0.3, 
which means that the column density increases with in¬ 
creasing orbital radius as S ex R 0 3 . The parameters are 
summarized in Table [4] 

Figure [H] shows the resultant polarization degree for 
the shallower temperature slope model. The calibration 
factor C of the formula CPqooj is set to be C = 1.3%. 
This figure shows that the formula CP 90 w can also fit 
the polarization fraction in the case of this model. P 9 quj 
is determined by grain size, and does not depend on the 
disk model. The all contribution of the change of the 
disk model comes into the calibration factor C. There¬ 
fore, we conclude that although the absolute value of 
the polarization fraction depends on the disk model, the 
constraints on the grain size do not so much depend on 
disk models; the polarization of the HL Tau can be re¬ 
produced only when the maximum grain size is around 
150 pm, although the minimum and maximum size of 
cimax can slightly depends on the disk model. 
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